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Abstract 

The production of J/y/ and </(2S) was studied with the ALICE detector in Pb-Pb collisions at the LHC. 
The measurement was performed at forward rapidity (2. 5 <y < 4) down to zero transverse momentum 
(Pt) in the dimuon decay channel. Inclusive J/y/ yields were extracted in different centrality classes 
and the centrality dependence of the average pi is presented. The i/xj/ suppression, quantified with 
the nuclear modification factor (Raa)^ was measured as a function of centrality, transverse momentum 
and rapidity. Comparisons with similar measurements at lower collision energy and theoretical models 
indicate that the J/\j/ production is the result of an interplay between color screening and recombination 
mechanisms in a deconfined partonic medium, or at its hadronization. Results on the </(2S) suppression 
are provided via the ratio of i/(2S) over J/ip measured in pp and Pb-Pb collisions. 
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1 Introduction 


At high temperature, lattice quantum chromodynamics predicts the existence of a deconfined phase of quarks 
and gluons where chiral symmetry is restored |[ll]. This state of matter is known as the Quark Gluon Plasma 
(QGP) 10], and its characterization is the goal of ultra-relativistic heavy-ion collision studies. 

Among the probes used to investigate the QGP and quantify its properties, quarkonium states are one of the 
most prominent and have generated a large amount of results both theoretical and experimental. According 
to the color-screening model i,li , measurement of the in-medium dissociation probability of the different 
quarkonium states could provide an estimate of the system temperature. Dissociation is expected to take 
place when the medium reaches or exceeds the critical temperature for the phase transition (Tc), depending 
on the binding energy of the quarkonium state. In the charmonium (cc) family, the strongly bound J/t/r 
could survive significantly above (1.5-2 T^) whereas Xc and t//^(2S) melting should occur near (1.1-1.2 
Tc) flSl. The determination of the in-medium quarkonium properties remains a challenging theoretical 
task. Intense and persistent investigations on the theory side are ongoing ||3]. Shortly after quarkonium 
suppression was suggested as a strong evidence of QGP formation, the first ideas of charmonium enhance¬ 
ment via recombination of c and c appeared lij^ Since then, the J/t/r enhancement mechanism has been 
more formalized and quantitative predictions ilOUldll were made. Since the charm quark density produced 
in hadronic collisions increases with energy il5ll . recombination mechanisms are predicted to give rise to a 
sizable J/y production at LHC energies, which is likely to partially compensate or exceed the J/t/r suppres¬ 
sion due to color-screening in the QGP. The observation of J/t// enhancement in nucleus-nucleus collisions 
via recombination would constitute an evidence for deconfinement and hence for QGP formation. In addi¬ 
tion, information for the characterization of the QGP can come from the study of the t/r(2S) meson, a state 
which is less strongly bound and not affected by higher mass charmonium decays with respect to the J/y- 
In the pure melting scenario, the relative production of V^(2S) with respect to J/y is expected to be very 
small at the LHC I^, which is not the case if recombination occurs ilhl lnl]. 


J/y suppression was observed experimentally in the most central heavy-nucleus collisions at the SPS lIlSL 


191], RHIC Il20l423n and LHC 12414281] . ranging from a center-of-mass energy per nucleon pair L/^ ) of 
about IVGeV to 2.76TeV. The t//^(2S) suppression was measured at the SPS ]29] and the LHC ]30]. The 
interpretation of these results is not straightforward as they are also subject to other effects, not all related to 
the presence of a QGP. A fraction of J/t/r originates from the strong and electromagnetic feed-down of the 
Xc and t//^(2S). Therefore, a melting of these higher mass states before they can decay into the J/y will lead 
to an effective suppression of the J/y yield already for a medium that does not reach the J/y dissociation 
temperature. Assuming charmonium states are initially produced with the same relative abundancies in 
Pb-Pb collisions as in pp collisions, the Xc and t/r(2S) melting would result in a reduction of the J/y yield 
of about 40% Oin . In addition, a non-prompt J/y and t//^(2S) component from the weak decay of beauty 
hadrons also contributes to the inclusive measurements. Since the beauty hadrons decay outside the QGP 
volume, this contribution is not sensitive to the color-screening of charmonia. Linally, a fraction of the 
J/y and t//^(2S) suppression can be ascribed to cold nuclear matter (CNM) effects, also present in proton- 
nucleus collisions ll32l. l33n. The CNM effects group together the nuclear absorption of the charmonia, the 
modification of the parton distribution functions (PDL) in the nuclei that leads to a reduction (shadowing) 
or an enhancement (anti-shadowing) of the cc pair production, and the energy loss of charm quarks in the 
nucleus. 

Numerous studies of J/y production in different collision systems at different energies are now available. 
Comparisons between experiments and to theoretical models can be made over wide kinematic ranges in 
rapidity and transverse momentum. We already published the centrality, transverse momentum (px) and 


rapidity (y) dependence of the J/y nuclear modification factor in Pb-Pb collisions at = 2.76TeV ||26L 
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2711 ■ In this paper, those results are extensively compared to available theoretical models and lower energy 
data. New results on the J/t/r (pj) and (pj) versus centrality, and on the centrality (pj) dependence of 
the J/y suppression for various pr (centrality) ranges are also presented. Furthermore, we show results on 
Vr(2S) in Pb-Pb collisions, measured via the [y{2S)/J/y] ratio, as a function of centrality. 


The remainder of this paper is organized as follows: the experimental apparatus and the data sample are 
presented in sections 2 and 3. Section 4 gives the definition of the observables used in the analysis. The 
analysis procedure is then described in sections 5 and 6. Systematic uncertainties are discussed in section 7. 
The J/y results are given in sections 8 and 9 while section 10 is dedicated to the '/r(2S) results. Finally, 
section 11 presents our conclusions. 


2 The ALICE detector 


The AFICE detector is described in detail in 13411 . At forward rapidity (2.5 <y < 4) the production of 
quarkonium states is studied in the muon spectrometer via their decay channels down to zero pj. 

In the AFICE reference frame, the positive z direction is along the counter-clockwise beam direction. The 
muon spectrometer covers a negative pseudo-rapidity (rj) range and consequently a negative y range. How¬ 
ever, due to the symmetry of the Pb-Pb system, the results are presented with a positive y notation, while 
keeping the negative sign for T]. 


The muon spectrometer consists of a ten-interaction-lengths (4.1 m) thick absorber, which filters the muons, 
in front of five tracking stations comprising two planes of cathode pad chambers each. The third station 
is located inside a dipole magnet with a 3 Tm field integral. The tracking apparatus is completed by 
a Muon Trigger system (MTR) composed of four planes of resistive plate chambers downstream from a 
seven-interaction-lengths (1.2 m) thick iron wall, which absorbs secondary hadrons escaping from the front 
absorber and low-momentum muons coming mainly from charged pion and kaon decays. A small-angle 
conical absorber protects the tracking and trigger chambers against secondary particles produced by the in¬ 
teraction of large rapidity primary particles with the beam pipe. Finally, a rear absorber protects the trigger 
chambers from the background generated by beam-gas interactions downstream from the spectrometer. 


In addition, the Silicon Pixel Detector (SPD), scintillator arrays (VO) and Zero Degree Calorimeters (ZDC) 
were used in this analysis. The SPD consists of two cylindrical layers covering |t]| < 2.0 and |t]| < 1.4 for 
the inner and outer ones, respectively, and provides the coordinates of the primary vertex of the collision. 
The VO counters, two arrays of 32 scintillator tiles each, are located on both sides of the nominal interaction 
point and cover 2.8 < rj <5.1 (VO-A) and —3.7 < rj < —1.7 (VO-C). The ZDC are located on either side 
of the interaction point atzRi±114m and detect spectator nucleons at zero degree with respect to the FHC 
beam axis. The VO and ZDC detectors provide triggering information and event characterization. 


3 Data sample 

The data sample analysed in this paper corresponds to Pb-Pb collisions at = 2.76TeV. These collisions 
were delivered by the FHC during 190 hours of stable beam operations spread over three weeks in November 
and December 2011. 

The Fevel-0 (FO) minimum bias (MB) trigger was defined as the coincidence of signals in VO-A and VO-C 
detectors synchronized with the passage of two crossing lead bunches. This choice for the MB condition 
provides a high triggering efficiency (> 95%) for hadronic interactions. To improve the trigger purity, a 
threshold on the energy deposited in the neutron ZDC rejects the contribution from electromagnetic dissoci¬ 
ation processes at the Fevel-1 (FI) trigger level. Beam induced background is further reduced at the offline 
level by timing cuts on the signals from the VO and the ZDC. 
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The charmonium analysis was carried out on a data sample, where in addition to the MB prerequisite, a 
trigger condition of at least one or two reconstructed muon candidate tracks in the MTR (trigger tracks) was 
required in each event. The MTR logic allows for programming several LO trigger decisions based on (i) 
the detection of one or two muon trigger tracks, (ii) the presence of opposite-sign or like-sign trigger track 
pairs and (iii) a lower threshold on the approximate transverse momentum (p*"®) of the muon candidates. 
The latter selection is performed by applying a cut on the maximum deviation of the trigger track from an 
infinite momentum track originating at the nominal interaction point. Due to the finite spatial resolution of 
the trigger chambers, this does not lead to a sharp cut in pj, and the corresponding p'"® threshold is defined 
in simulation as fhe px value for which fhe muon frigger probabilify is 50%. The following muon-specific 
LO friggers were used: 

- Single muon low pj (p*"® = 1 GeV/c): MSL 

- Opposife-sign dimuon low pj (p*"® = 1 GeV/c on each muon): MUL 

- Like-sign dim non low pj (px'® = IGeV/c on each muon): MLL 


A dafa sample of 17.3 • 10^ Pb-Pb collisions was collecfed wifh fhe pp-MB frigger, defined as fhe coinci¬ 
dence of fhe MB and MUL conditions. A scaling factor Lnorm is compufed for each run — corresponding 
fo a few hours maximum of continuous dafa faking — in order to normalize fhe number of pp-MB friggers 
to fhe number of equivalenf MB friggers. If is defined as fhe ratio, in a MB dafa sample, befween fhe fofal 
number of evenfs and fhe number of evenfs fulfilling fhe pp-MB frigger condition. If should be nofed fhaf 
fhe MB sample used in fhis calculation was recorded in parallel fo fhe pp-MB friggers. The Lnorm value, 
30.56 ± 0.01 (sfaf.) ± 1.10(sysl.), is given by fhe average over all runs weighted by fhe sfafisfical uncer- 
fainfies. A small fraclion of opposite-sign dimuons were misidenfified by fhe frigger algorifhm as like-sign 
pairs. Allhough for fhe J/tp if amounls fo less fhan 1% when considering fhe full sample, if increases up 
to 4% af high px in peripheral collisions. In fhis analysis, fhe missing fraction of opposife-sign dimuons 
was recovered by exlracling fhe number of produced J/tp and tp(2S) from fhe union of fhe MUL and MLL 
dafa sample (MULUMLL). This is differenf from fhe selection applied in fhe former paper 12711 . where only 
fhe MUL dafa sample was used. On fhe olher hand, fhe efficiency of fhe frigger algorifhm fo delermine 
fhe sign of fhe muon pairs does nol impacf fhe normalizalion of fhe collected dafa sample fo fhe number of 
equivalenf MB evenfs described above. This was cross-checked by computing fhe normalization facfor of 
fhe MULUMLL dafa sample, resulting in less fhan 1 % difference in fhe exlracled number of equivalenf MB 
evenfs. 


The infegrafed luminosity corresponding fo fhe analysed dafa sample is • Lnorm/<7pb-Pb = 

68.8±0.9(sfaf.) ±2.5(sysf. tyiorm)^ 45 (sysf. apb-pb) qb^' using an inelastic Pb-Pb cross section apb-pb = 
7.7±0.llJ];^b il. 


4 Definition of observables 


The cenfralify determination is based on a fit of the VO signal amplitude distribution as described in 13611 . 
Variables characterizing the collision such as the average number of participant nucleons ((Apart)) and the 
average nuclear overlap function ((7aa)) for each centrality class are given in Tab. [T] In this analysis a cut 
corresponding to the most central 90% of the inelastic nuclear cross section was applied as for these events 
the MB trigger is fully efficient and the residual contamination from electromagnetic processes is negligible. 
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Centrality 

(Apart) 

{Taa) (mb^i) 

Centrality 

(Apart) 

{Taa) (mb^*) 

0 - 10 % 

356.0±3.6 

23.44±0.76 

0 - 20 % 

308.1±3.7 

18.91±0.61 

10 - 20 % 

260.1±3.8 

14.39±0.45 

0-40% 

232.6±3.4 

12.88±0.42 

20-30% 

185.8±3.3 

8.70±0.27 

0-90% 

124.4±2.2 

6.27±0.21 

30-40% 

128.5±2.9 

5.00E0.18 

20 ^ 0 % 

157.2±3.1 

6.85±0.23 

40-50% 

84.7±2.4 

2 . 68 ± 0.12 

20-60% 

112 . 8 ± 2.6 

4.42±0.16 

50-60% 

52.4E1.6 

1.317±0.071 

40-60% 

68 . 6 ± 2.0 

1.996±0.097 

60-70% 

29.77±0.98 

0.591±0.036 

40-90% 

37.9±1.2 

0.985±0.051 

70-80% 

15.27±0.55 

0.243±0.016 

50-90% 

26.23±0.84 

0.563±0.033 

80-90% 

7.49±0.22 

0.0983±0.0076 

60-90% 

17.51±0.59 

0.311 ± 0.020 


Table 1: The average number of participant nucleons (Apart) and the average value of the nuclear overlap function 
(Taa) with their associated systematic uncertainties for the centrality classes, expressed in percentages of the nuclear 
cross section 1^ . used in these analyses. 


For each centrality class i, the measured number of J /y is normalized to the equivalent number of 

minimum bias events To obtain one simply multiplies the number of /i/r-MB triggered 

events by the Enorm factor scaled by the width of the centrality class. Corrections for the branching ratio of 
the dimuon decay channel = 5.93 ±0.06%) and for the acceptance times efficiency (A x e') 

of the detector are then applied. The J/y yield (Tj^^) in a centrality class i is given by 


j2yi 

dpxdy 




^j/vz/d/^Tdy 


'-^events ^ 


( 1 ) 


It is then combined with the inclusive J/t/r cross section measured in pp collisions at the same energy to 
form the nuclear modification factor Raa defined as 


f^AA(PTA) = 


d^^j/y/dPrdy 

(rAA)'-d2ajJ^/dpTdy' 


( 2 ) 


The Pt and y integrated J/t/r cross section is < 8 GeV/c, 2.5 < y < 4) = 3.34 ± 0.13(stat.) ± 

0.24(syst.) ±0.12(luminosity)lj 07 (polarization)p,b s. 

The ALICE measurements reported here refer to inclusive J/t/r yields, i.e. include prompt J/t/r (direct J/yr 
and feed-down from t//^(2S) and Xc) and non-prompt J/t/r (decay of B-mesons). Contrary to prompt J/t/r, 
J/t/r from B-meson decays do not directly probe the hot and dense medium created in the Pb-Pb collisions. 
Beauty hadron decays occur outside the QGP, so the non-prompt J/y/ Raa is instead related to the energy 
loss of the beauty quarks in the medium. Although the prompt i/Y Raa cannot be directly measured with 
the ALICE muon spectrometer, it can be evaluated via 


R 


prompt 

AA 


Ty 77 onon-prompt 

^AA • ^aA 


(3) 


where Eb is the fraction of non-prompt to inclusive J/t/r measured in pp collisions, and prompt 

clear modification factor of J/t/r from B-meson decays in Pb-Pb collisions. The non-prompt and prompt 

J/t/r differential cross sections as a function of and y were measured by LHCb in pp collisions at 
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y/s =2.76 and 7 TeV OSL 13911 in a kinematic range overlapping with that of the ALICE muon spectrom¬ 
eter. Therefore, one can extract the and y dependence of Fq from these data and use it in Eq. [3] A 

presents further complications. We have thus chosen two extreme hy- 


reliable determination of R 


,non-prompt 
AA 


potheses, independent of centrality, corresponding to the absence of medium effects on beauty hadrons 
= 1) or to a complete suppression (7 ?a*a = 0), to evaluate conservative limits on 

An excess of J/t// compared to the yield expected assuming a smooth evolution of the J/t/r hadro-production 
and nuclear modification factor was observed in peripheral Pb-Pb collisions at very low pj i40n . This excess 
might originate from the photo-production of J/t/r. This contribution is negligible in pp collisions — from 
LHCb measurement at y/s = 7TeV |^], it is ^(10^^)% — but it is enhanced by a factor ^(10“^) in Pb-Pb 
collisions, thus reaching the order of magnitude of the observed excess. The J/y/ coherent photo-production 
has been measured in ultra-peripheral Pb-Pb collisions ll^ . It is centered at very low pj, with ~ 98% of 
these J/y/ below 0.3 GeV/c. An incoherent photo-production component is also observed in ultra-peripheral 
Pb-Pb collisions. About 30% of this contribution has a pj < 0.3 GeV/c, the rest being mainly located in 
the Pj range 0.3-1 GeV/c. The influence of possible photo-production mechanisms on the inclusive J/yr 
Raa presented in this paper has been evaluated by repeating the analysis placing a low pj threshold on the 
J/yr of 0.3 GeV/c. Assuming that the observed excess in peripheral Pb-Pb collisions is indeed due to the 
photo-production of J/yr, and that the relative contribution of the incoherent over coherent components is 
the same as the one estimated in ultra-peripheral collisions, then this selection would remove about 75% of 
the full photo-production contribution. Numerical values of Raa with the low pj threshold at 0.3 GeV/c are 
given in the Appendix All the figures and values presented in the paper refer to the inclusive J/yr Raa but 
estimates of the difference between the inclusive and hadronic (without J/yr photo-production) J/yr Raa, 
are indicated where appropriate. 

The results for the y/^(2S) analysis are given in terms of the ratio of their production cross sections (or, 
equivalently, of their production yields), expressed as 


y/(2S)/J/yr = 


^W{2S) (A ^ 


N\/^ (Axe%( 2 s)' 


(4) 


When forming such a ratio the normalization factor Ngvents cancels out, as do most of the systematic uncer¬ 
tainties on A X e corrections. The double ratio [yr(2S)/J/yr]pj^p(, / [yr(2S)/J/yr]pp is used in order to directly 
compare the relative abundances of yr(2S) and J/yr in nucleus-nucleus and pp collisions. 


5 Signal extraction 

After a description of the muon selection procedure, we present here the two methods used to extract the 
J/yr and y/^(2S) signals. The first one is directly based on fits of the invariant mass distribution while 

the second one makes use of the event mixing technique to subtract the combinatorial background. 

5.1 Muon reconstruction 

The muon reconstruction starts with the exclusion of parts of the detector that show problems during data 
taking such as high voltage trips, large electronic noise, pedestal determination issues. This selection is 
performed on a run-by-run basis to account for the time evolution of the apparatus. After pedestal sub¬ 
traction, the adjacent well-functioning pads of both cathodes of each tracking chamber having collected a 
charge are grouped to form pre-clusters. These pre-clusters might be the superposition of several clusters 
of charges deposited by several particles crossing the detector close to each others. The number of clusters 
of charges contributing to the pre-cluster and their approximate location are determined with a Maximum 
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Likelihood - Expectation Maximization (MLEM) algorithm. It assumes that the charge distribution of each 
single cluster follows a two-dimensional integral of the Mathieson function 1430. If the estimated number 
of clusters is larger than 3, the pre-cluster is split into several groups of 1, 2 or 3 clusters selected with the 
minimum total coupling to all the other clusters into the pre-cluster. Each group of clusters is then htted 
using a sum of Mathieson functions, taking the MLEM results as a seed, to extract the precise location of 
where the particles crossed the detector. The overall spatial resolution is around 200 (550) p.m in average in 
the (non-)bending direction. 


The track reconstruction starts from the most downstream stations, where the multiplicity of secondary 
particles is smallest, by forming pairs of clusters in the two chambers of station 5(4), and deriving the pa¬ 
rameters and associated errors of the resulting muon track candidates. The candidates are then extrapolated 
to the station 4(5), validated if at least one compatible cluster is found in the station and duplicate tracks 
are removed. The procedure continues extrapolating the tracks to stations 3, 2 and 1, validating them by 
the inclusion of at least one cluster per station. The selection of compatible clusters is based on a 5a cut 
on a computed from the cluster and track local positions and errors. If several compatible clusters are 
found in the same chamber, the track is duplicated to consider all the possibilities and for each of them 
the track parameters and associated errors are recomputed using a Kalman hlter. At each of the tracking 
steps, the track candidates, whose parameters indicate that they will exit the geometrical acceptance of the 
spectrometer in the next steps are removed. At the end of the procedure, the quality of the track is improved 
by adding/removing clusters based on a 4a cut on the local x^ and fake tracks sharing clusters with others 
in the three outermost stations with respect to the interaction point are removed. The choice of the x^ cuts 
is a compromise between maximizing the tracking efficiency (< 1-2% muon rejection) and minimizing the 
amount of fake tracks (negligible background for this analysis). Einally, muon track candidates are extrap¬ 
olated to the interaction vertex measured by the SPD taking into account the energy loss and the multiple 
Coulomb scattering in the front absorber. 


An accurate measurement of the tracking chamber alignment is essential to reconstruct the tracks with 
enough precision to identify resonances in the invariant mass spectrum, especially the t//^(2S) for 

which the signal-to-background ratio is low. The absolute position of the chambers was hrst measured 
using photogrammetry before the data taking. Their relative position was then precisely determined using 
a modihed version of the MILLEPEDE package 14411 . combining several samples of tracks taken with and 
without magnetic held. The small displacement of the chambers when switching on the dipole was measured 
by the Geometry Monitoring System (an array of optical sensors hxed on the chambers) and taken into 
account. The resulting alignment precision is ~ 100 pm, leading to a reconstructed J/y/ invariant mass 
resolution of about 70 MeV/c^, and about 10% higher for the V^(2S). The resolution is dominated by the 
energy loss huctuation and multiple Coulomb scattering of the muons in the front absorber. More details on 
the muon spectrometer performances are given in 14511 . 


In this analysis, the muon track candidates also have to fulhll the following requirements. Eirst, the recon¬ 
structed track must match a trigger track with a above the threshold set in the MTR for triggering the 
event (IGeV/c in this analysis). The trigger track is reconstructed from the average position of the hred 
strips on the two trigger stations, as computed by the trigger algorithm. The matching is based on a 4a cut 
on a x^ computed from the tracker and trigger track parameters and errors including the angular dispersion 
due to the multiple Coulomb scattering of the muon in the iron wall. Second, the transverse radius coordi¬ 
nate of the track at the end of the front absorber must be in the range 17.6 < Rahs < 89.5 cm. Muons exiting 
the absorber at small and large angles, thus outside the Rahs cut range, have crossed a different amount of 
material, either the beam shield or the envelope of the absorber, affecting the precision of the energy loss 
and multiple Coulomb scattering corrections. Third, in order to remove muon candidates close to the edge 
of the spectrometer acceptance, a cut on the track pseudo-rapidity —4 < Tj < —2.5 is applied. 
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5.2 J/t/r signal 

J/Y candidates are formed by combining pairs of opposite-sign tracks reconstructed within the geometrical 
acceptance of the muon spectrometer. The aforementioned cuts at the single muon track level remove most 
of the hadrons escaping from or produced in the front absorber, as well as a large fraction of low pj muons 
from pion and kaon decays, secondary muons produced in the front absorber, and fake tracks. The J/y peak 
becomes visible in the invariant mass spectra even before any background subtraction. At the dimuon 

level only cuts on rapidity (2.5 < y < 4) and transverse momentum {pj < 8 GeV/c) are applied. The J/y 
raw yields are extracted by using two different methods. 


In the first method, the opposite-sign dimuon invariant mass distribution is fitted with a sum of two functions. 
The signal is describe d by a double-sided Crystal Ball function (CB2). This function is an extension of the 
Crystal Ball function 14611 . i.e. a Gaussian with a power-law tail in the low mass range, with an additional 
independent power-law tail in the high mass range. The CB2 function reproduces very well the J/y line 
shape in the Monte Carlo (MC) simulations. The underlying continuum is described by a variable width 
Gaussian function. This function is built on a Gaussian form, whose width is dependent on the invariant 
mass of the dimuon. It was checked that including or excluding a y{'^^) contribution in the fitting procedure 
has a negligible effect on the extracted J/y yield within the present statistical and signal-extraction-related 
systematic uncertainties. Since the significance of the y{^^) signal in the centrality, pj and y intervals used 
for the J/y analysis is too small to extract its contribution, we do not include it in the fit for this analysis. 
During the fitting procedure, the width of the J/t/r peak is kept as a free parameter as it cannot be reproduced 
perfectly in simulations, and its value varies from 65 to 76 MeV/c^ (larger than those from MC by about 
5-10%). The pole mass is also kept free although the differences observed between data and simulation are 
at the per mille level. The tail parameters cannot be constrained by the fit. Therefore they are fixed fo values 
obtained from an embedding simulation (described in section and adjusted for each pj and y interval 
under study in order to account for the observed dependence on the J/y kinematics. On the contrary, the 
J/y shape does not show a dependence on centrality, hence the CB2 tail parameters tuned on a centrality 
integrated MC sample are used in all the bins. Figure [T] presents fits of the opposite-sign dimuon invariant 
mass (m^^) distributions for different pr ranges in central (top row) and peripheral (bottom row) collisions. 
The signal-to-background ratio (S/B) and the significance (S/\/S -|-B) of the signal are evaluated within 3 
standard deviations with respect to the J/y pole mass. The S/B varies from 0.2 to 6.5 when going from 
the most central collisions to the most peripheral ones. Integrated over centrality and y (pr), the S/B ranges 
from 0.2 (0.2) to 1.2 (0.6) with increasing px (y)- In all the centrality, px or y intervals considered in this 
analysis, the significance is always larger fhan 8. 


In the second method, the combinatorial background is subtracted using an event-mixing technique. The 
opposite-sign muon pairs from mixed-events are formed by combining muons from single muon low px 
(MSL) triggered events. In order to limit the effect of efficiency fluctuations between runs and to take into 
account the dependence of muon multiplicity and kinematic distributions on the collision centrality, events 
in the same run and in the same centrality class are mixed together. The mixed-event spectra are normalized 
to the data using the combination of the measured like-sign pairs such as 




dm 


MM 


= 


’ /dA++dA__ , 

^\l dm,, 


‘MM 


(5) 


where A^+-, N++ and N _are the number of opposite-sign, positive like-sign and negative like-sign muon 
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Fig. 1 : Fit to the opposite-sign dimuon invariant mass distribution in the 0-20% (upper row) and 40-90% (lower row) 
centrality classes, for 2.5 <y < 4, in various pj intervals. 


pairs. The R factor in Eq.[5]is defined by 


R = 


djym^xed 

dnififi 


2 


dN'T 


dmu 


^mixcd 

dnififi 


( 6 ) 


and is introduced in order to correct for differences in acceptance between like-sign and opposite-sign muon 
pairs. Above a dimuon invariant mass of 1.8 GeV/c^, the R factor is equal to unity with deviations smaller 
than 1%. The accuracy of the mixed-event technique was assessed by comparing the distributions of like- 
sign muon pairs from mixed-events to the same-event ones, which agree within 1% over the mass, and 
y ranges under study. This agreement justifies the use of the normalization given by Eq. |5l which implies 
that the correlated signal in the like-sign dimuon spectra is negligible with respect to the combinatorial 
background. The mass spectra of the opposite-sign mixed-event pairs are then subtracted from the data. 
The resulting background-subtracted spectra are fitted following the same procedure as in the first method, 
except that the variable width Gaussian function is replaced by an exponential function to account for resid¬ 
ual background. Eigure [2] shows fits of the background-subtracted opposite-sign dimuon invariant mass 
distributions for different pr ranges in central (top row) and peripheral (bottom row) collisions. 


5.3 t/^(2S) signal 

The invariant mass spectra used to extract the [t//^(2S) /J/y] rati® obtained in the same way as described in 
the previous section, implementing the same cuts applied at the muon and dimuon levels. In order to improve 
the significance of the t//^(2S) signal, a wider centrality and pj binning than the one used for the J/y analysis 
was adopted, and the analysis is performed integrated over the full rapidity domain 2.5 < y < 4. The fits 
to the invariant mass spectra are performed by modeling the t//^(2S) signal with a CB2 function. Given the 
very low S/B ratio, the normalization is chosen as the only free parameter for vr(2S). The tails of the CB2 
function describing the t//^(2S) are fixed fo fhose exfracfed for fhe J/t/r. The posifion of fhe V^(2S) pole 
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Fig. 2: Fit to the opposite-sign dimuon invariant mass disttibution after background subtraction in the 0-20% (upper 
row) and 40-90% (lower row) centrality classes, for 2.5 <y < 4, in various intervals. 


mass is fixed to the one of the J/y/, shifted by the corresponding Am = value taken from the 

PDG 114711 . The width of the y/^(2S) is fixed fo fhe one of fhe J/yr scaled by fhe rafio <7||/(2 S)/<^j/v/ estimated 
from MC simulafions. 


Fils of fhe invarianl mass speclra showing fhe yr(2S) are visible in Fig.[3]for fhe pj < 3GeV/c interval in 
the centrality classes 20^0%, 40-60% and 60-90%. For the other intervals in centrality and pj the y/^(2S) 
signal could not be extracted, i.e. the ratio [y/^(2S)/J/y^] is consistent with zero. In these cases, only the 95% 
confidence level upper limil is computed. 


6 Acceptance and efficiency correction 


In fhe J/yr analysis, embedding simulafions are used fo compule fhe cenlralily, pj and y dependences of 
the acceptance times efficiency factor (A x e). The Monfe Carlo embedding lechnique consisls of adding 
the detector response from a simulated signal event (charmonium in this case) to a real data event, and then 
performing the reconstruction as for real events. This has the advantage of providing the most realistic 
background conditions, which is necessary for Pb-Pb collisions where high multiplicities are reached: at 
rj = 3.25, dAch/drj 1450 for the 0-5% most central events 14811 . This leads to a large detector occupancy, 
which can reach about 3% in the most central collisions and alter the track reconstruction efficiency. 

Monfe Carlo J/yr were embedded in MB friggered evenfs recorded in parallel fo fhe opposite-sign dim non 
friggered evenfs. Only one J/yr was simulated per even! af fhe position of fhe real even! primary vertex 
reconsfrucfed by fhe SPD. The shapes of fhe inpuf MC px and y disfribufions were funed fo mafch fhe 
measured disfribufion in Pb-Pb collisions (see discussion in section \T72\ . The muons from fhe J/y/ de¬ 
cay were fhen fransported fhrough a simulation of fhe AFICE detector using GEANT3 14911 . The defecfor 
simulafed response was fhen merged wifh fhaf of a real Pb-Pb even! and fhe resulf was processed by fhe 
normal reconsfrucfion chain. Embedding simulafions were performed on a run-by-run basis fo accounf for 
fhe fime-dependenf sfafus of fhe fracking chambers. The residual misalignmenf of fhe defecfion elemenfs. 
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Fig. 3: Opposite-sign dimuon invariant mass distribution for the 20^0%, 40-60% and 60-90% centrality classes, for 
2.5 < y < 4 and pj <3 GeV/c, before (top row) and after background subtraction (bottom row) via event mixing. In 
these intervals the i/(2S) signal is extracted whereas in all other centrality and pj intervals, only the 95% confidence 
level upper limits are provided. 


whose amplitude is evaluated by analyzing the residual distance between the clusters and the tracks in data, 
was also taken into account. For the trigger chambers, the efficiency maps measured in data were used in 
the simulations. 

The left panel of Fig.|4]shows the J/t// A x e as a function of collision centrality in the rapidity domain 2.5 < 
y < 4 and in the pj range pj < 8GeV/c. We observe a relative decrease of 8% of the J/t/r reconstruction 




Fig. 4: The J/y/ acceptance times efficiency, shown as a function of centrality (left) and as a function of pj and y for 
the centrality class 0-90% (right). The vertical error bars in the left panel represent the statistical uncertainties. 
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efficiency from the 80-90% centrality class to the 0-10% centrality class. This decrease is mostly due to a 
drop of about 3% of the single muon trigger efficiency in the most central collisions whereas the decrease of 
the single muon tracking efficiency is only on the order of 1%. When considering specific pj or y intervals, 
a maximum relative variation of ~ 30% of the A x e decrease with centrality is observed. The right panel of 
Fig. m shows the pj versus y dependence of A x e. The rapidity dependence of A x £ reflects the geometrical 
acceptance of the muon pairs with a maximum centered at the middle of the rapidity interval and a decrease 
towards the edges of the acceptance. The px dependence of A x £ is non-monotonic, with a minimum at 
Pt ~ 1-8 GeV/c corresponding to J/i/r kinematics for which one of the decay muons does not fall into the 
muon spectrometer acceptance. 

For the V^(2S) resonance, the embedding technique was not used. Since, in this case, only the ratio 
[y{2S)/J/y] is extracted, the Ax£ correction factors for both resonances were evaluated through pure 
signal MC simulations, assuming that the dependence of the efficiency as a function of the centrality is the 
same for J/i/r and i/r(2S), and therefore cancels out in the ratio. The effect of possible differences in the 
centrality dependence of A x £ was studied and included as a source of systematic uncertainty. 

7 Systematic uncertainties 

In the following, each source of systematic uncertainty is detailed. Most of them affect the J/yr and V^(2S) 
results identically and vanish in the [i/r(2S)/J/vr] ratio. Systematic uncertainties specific to the '/^(2S) 
analysis are explicitly mentioned. 

7.1 Signal extraction 

The systematic uncertainty on the signal extraction results from several fits of the invariant mass spectra, 
where signal line shape parameters, background description and fit range are varied as detailed below. In 
each centrality, px and y intervals, the raw yield and the statistical uncertainty are given by the average of 
the results obtained from the different fits. The corresponding systematic uncertainty is defined as the RMS 
of these results. It was also checked that every individual result differs from the mean value by less than 
three RMS. 

The J/y line shape is well described by the CB2 function, whose pole mass and width are constrained 
by the data while the tail parameters have to be fixed to values extracted from the embedding simulation. 
Alternatively, another set of tails was extracted from pp data, where a large statistics and a better S/B were 
available. In this case, the px and y dependence of the tail parameters could not be determined with sufficient 
precision, so the same values were used for all px and y intervals. In the event mixing approach, the 
influence of different normalizations of the opposite-sign mixed-event spectrum to the opposite-sign same- 
event spectrum was investigated. We have tested a normalization performed on a run-by-run basis or after 
merging of all the runs, and a normalization based on the integral of the invariant mass spectrum in the 
intermediate mass region (1.5 < < 2.5 GeV/c^). None of these tests showed deviations larger than 

1% in the number of extracted J/i/r, and thus were not included in the tests used to extract the systematic 
uncertainty on the signal extraction. The fit range of the invariant mass spectra was also varied considering a 
narrow (2.3 < < 4.7GeV/c^) and a wide (2 < < 5GeV/c^) interval. Finally, all the combinations 

of signal line shape, background description (with or without using the event-mixing technique) and fit range 
are performed to account for possible correlations. 

The same procedure as above was applied when the V^(2S) signal was included in the fit function for the 
specific centrality and px intervals presented in this analysis. To account for the fact that the V^(2S) width 
was fixed to the one of the J/y scaled by the ratio <7|^(2S)estimated from MC simulations, all the fits 
were repeated varying the scaling factor by ±10%. This variation accounts for the fluctuations observed in 
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pp data when fitting the invariant mass spectra leaving the width of the t//^(2S) free or fixing if as described 
above. 

The sysfemafic uncerfainfy on fhe signal exfracfion varies wifhin fhe 1^% range depending on fhe cenfralify 
class. Considering fhe pj infervals 0-2, 2-5 and 5-8GeV/c used for fhe Raa mulli-differenlial sfudies, we 
obfain sysfemafic uncertainties in fhe ranges 1^%, 1-4% and 1-3%, respectively. As a funcfion of p^, the 
systematic uncertainty on the signal extraction varies from 1% to 4%; for the centrality 0-20%, 20-40% 
and 40-90%, the values are in the ranges 1-5%, 1^% and 1-2%, respectively. As a function of y, the 
systematic uncertainty on the signal extraction varies from 1% to 4%. Concerning the t//^(2S) analysis, in 
the intervals where the signal was extracted, the systematic uncertainty is 14%, 45% and 24% for centrality 
ranges 60-90%, 40-60% and 20-40% for pj < 3GeV/c. 

7.2 Monte Carlo input parametrization 

The estimation of A x e factors depends on the charmonium pj and y shapes used as input distributions in the 
MC simulation. In order to evaluate the sensitivity of the results on this choice, several MC simulations were 
performed, each one including modified pj and y disfribufions. For fhe J/t/r, the modification of the shapes 
was done in order to take into account the possible correlation between p^ and y (as observed by LHCb in 
pp collisions and the correlation between p-^ (y) and the centrality of the collision (as reported in this 
paper). A systematic uncertainty of 3% is found for A x e integrated over p^ and y and is taken as correlated 
as a function of the centrality. The pj (y) dependence of this uncertainty varies in the range 0-1% (3-8%). 
The larger effect seen in the y dependence occurs at the low and high limits, where the acceptance falls 
steeply. 

The same procedure was followed for the V^(2S), assuming that the correlations between pj and y and with 
the centrality are of the same magnitude as those observed for the J/t/r. A systematic uncertainty of 2% is 
evaluated for the [v^(2S) /^/y\ in th® Pi < 3GeV/c interval. 

7.3 Centrality dependence of the [t/^(2S)/J/v^] Axe 

The embedding technique was not used for the V^(2S) MC simulations as we have assumed the same Axe 
dependence as a function of the centrality for the t//^(2S) and the J/t/r. In order to evaluate the systematic 
uncertainty introduced by this assumption, a conservative ±30% variation of the Axe loss as a function of 
centrality was applied to the t//^(2S). This corresponds to the maximum variation of the Axe loss between 
peripheral and central collisions observed for the J/t/r in different pj and y intervals. The effect on the 
(AX e)j/^/ (Ax £)y( 2 s) is 1% or lower in all the centrality classes considered. Since this effect is much 
smaller than the systematic uncertainty on the signal extraction, it is neglected. 

7.4 Tracking efficiency 

The tracking algorithm, as described in section [5Al does not require all the chambers to have fired to recon¬ 
struct a track. This redundancy of the tracking chambers can be used to measure their individual efficiencies 
from data, and since such efficiencies are independent from each other, we can combine them to assess the 
overall tracking efficiency. This evaluation of the tracking efficiency is not precise enough to be used to 
directly correct the data, because only the mean efficiency per chamber can be computed with the statistics 
available in each run. However, by comparing the result obtained from data with the same measurement 
performed in simulations, we can control the accuracy of these simulations and assess the corresponding 
systematic uncertainty on the Axe corrections. 

A 9% relative systematic uncertainty is obtained for the J/v^ by comparing the measured tracking efficiency 
in simulations and in peripheral Pb-Pb collisions. This uncertainty is constant and fully correlated as a 
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function of centrality. From low to high px Cy)^ the systematic uncertainty varies from 9% to 7% (7% to 
6% with a maximum of 12% at y ~ 3.25). On top of that, a small difference was observed in the centrality 
dependence of this measurement between data and embedding simulations. This results in an additional 1 % 
systematic uncertainty in the 0-10% centrality class and 0.5% in 10-20%. 

Another systematic uncertainty can arise from correlated dead areas located in front of each other in the same 
station, which cannot be detected with the method detailed above. A dedicated study has shown that this 
effect introduces a 2% systematic uncertainty, fully correlated as a function of centrality and predominantly 
uncorrelated as a function of pj and y. 

In the [t/r(2S)/J/vr] ratio the systematic uncertainty on the tracking efficiency largely cancels out because 
the V^(2S) and J/t/r decay muons have similar pj and y distributions and, therefore, cross about the same 
regions of the detector. Since the possible remaining systematic uncertainty is much smaller than that on the 
signal extraction, it is neglected in this analysis. 


7.5 Trigger efficiency 

The systematic uncertainty on the J/t/r Axe corrections related to the trigger efficiency has fwo origins: 
the intrinsic efficiency of the trigger chambers and the response of the trigger algorithm. The first part was 
determined from the uncertainties on the trigger chamber efficiencies measured from dafa and applied fo 
simulafions. Propagating fhese efficiencies in J/t/r simulations resulfs in a 2% sysfemafic uncerfainfy on fhe 
Axe correcfions, fully correlafed as afuncfion of cenfralify and mainly uncorrelafed as afuncfion of pj and 
y. The effecf of fhe sysfemafic uncerfainfy on fhe shape of fhe trigger response as a function of the muon px 
was determined by weighting MC J/y/ decay muons with different trigger response functions obtained in 
data and simulations. These functions were defined as fhe fraction, versus px, of fhe single muons passing 
a 0.5GeV/c Px”® fhreshold fhaf also satisfy fhe IGeV/c Px”® fhreshold used in fhis analysis. The resulting 
sysfemafic uncerfainfy on fhe J/t/r Axe correction infegrafed over px and y is 1%. As a funcfion of px, 
if amounfs fo 3% for px < IGeV/c and 1% elsewhere. As a funcfion of y, a 1% uncorrelafed sysfemafic 
uncerfainfy was obfained. These uncerfainfies are fully correlafed as a funcfion of cenfralify. 

The sysfemafic uncerfainfy on fhe modificafion of fhe frigger response as a funcfion of cenfralify, i.e. for 
increasing mulfiplicify, was assessed by changing fhe defecfor response (space size of fhe deposifed charge) 
fo fhe passage of particles in embedding simulations. The corresponding uncerfainfies on the J/t/r Ax£ 
corrections are 1% in the 0-10% and 10-20% centrality classes, and 0.5% in 20-30% and 30-40%. 

As for the case of tracking efficiency, fhis source of sysfemafic uncerfainfy largely cancels ouf in fhe 
[t/r(2S) /J/y] ratio is neglecfed. 

7.6 Matching efficiency 

The sysfemafic uncerfainfy on fhe mafching efficiency befween the tracking and the trigger tracks is 1%. It 
is given by the differences observed between data and simulations when applying different cuts on the 
matching between the track reconstructed in the tracking chambers and the one reconstructed in the trigger 
chambers. This uncertainty is fully correlated as a function of the centrality and largely uncorrelated as a 
function of px and y. 

Also in this case, the effect on the [t/r(2S)/J/t/r] ratio is negligible. 


7.7 pp reference 


The statistical and systematic uncertainties on the measurement of the J/y differential cross section in 
pp collisions at ^/s = 2.76TeV are available in 0711 . The statistical uncertainty is combined with that of 
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the Pb-Pb measurement when calculating the /?aa as a function of pj and y, but is considered as a fully 
correlated systematic uncertainty as a function of the centrality. The correlated and uncorrelated part of the 
systematic uncertainty on the pp reference as a function of pj and y are both fully correlated as a function 
of the centrality. 

The V^(2S) statistics in the y/s = 2.76TeV pp data sample are too low to be used for the normalization 
of the [vr(2S)/J/t/r]pj^p(^ ratio. For this reason, pp results obtained at higher energy {y/s = 7TeV) if^ 
were used, thus introducing an additional source of systematic uncertainty. An interpolation procedure, as 
the one described in iBill . was applied in order to extract the [t/r(2S)/J/t/r]pp ratio at y/s = 2.76TeV. The 
discrepancy between the result of this interpolation in the kinematic range pr < 3GeV/c 2.5 < y < 4 and 
the value obtained at y/s = 7 TeV is 10%: this relative difference is included in the systematic uncertainty 
on the pp reference. 

7.8 Normalization 

The systematic uncertainty on the normalization is the one attached to the scaling factor Enonn and amounts 
to 4%. This value corresponds to one standard deviation of the distribution of the Enorm computed for each 
run used in the analysis. This systematic uncertainty is fully correlated as a function of the centrality, pj and 
J- 


7.9 Others 


Systematic uncertainties on the nuclear overlap function (Eaa) are available in Tab. [T] Another systematic 
uncertainty on the definition of the centrality classes arises from the VO amplitude cut, which corresponds 
to 90% of the hadronic cross section li^ . A maximum uncertainty of 5% is obtained in the centrality class 
(80-90%) vanishing with increasing centrality or in wider centrality classes. 


Systematic uncertainties due to the unknown polarization of the J/t/r are not propagated and we assume 
that J/y production is unpolarized both in pp and in Pb-Pb collisions. In pp collisions at y/s = 7 TeV, 
J/Y polarization measurements at mid-rapidity (px > lOGeV/c) and forward-rapidity (px > 2GeV/c) are 
compatible with zero 15114531] . In Pb-Pb collisions, J/t/r mesons produced from initial parton-parton hard 
scattering are expected to have the same polarization as in pp collisions and those produced from charm 
quarks recombination in the medium are expected to be unpolarized. 


7.10 Summary 

The systematic uncertainties related to the J/y analysis are summarized in Tab. |2] Concerning the t//^(2S) 
analysis, most of the systematic uncertainties cancel out in the [v/^(2S)/J/t/r] ratio and the main contributors 
are the signal extraction (14^5%) and the pp reference (10%). 


8 Inclusive J/y mean transverse momentum 

The px dependence of the J/y yields per MB collision, defined by Eq. [H was sfudied for fhree cenfralify 
classes (0-20%, 20-40% and 40-90%) and is displayed in Fig. [5] The sfafisfical uncerfainfies appear as 
vertical lines. The sysfemafic uncertainties uncorrelafed as a funclion of px are shown as open boxes, while 
fhe ones fully correlafed as a function of px buf uncorrelafed as a function of centralify are shown as shaded 
areas (mosfly hidden by fhe poinfs). The global sysfemafic uncertainfy, fully correlafed as a funcfion of 
cenfralify and px, is quofed direcfly in fhe figure. Numerical values for fhe J/y yields can be found in 
Appendix 13 The inclusive J/y mean fransverse momenfum was compufed by filling fhe px dislribulion of 
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Sources 

Centrality 

Pt < 8 GeV /c [271 p^ bins 

0-90% [27J 

Pt 

centrality bins 

y[22J 

Signal extraction 

1-3 

1-4 

1^ 

1-5 

1^ 

MC parametrization 

3* 

1-3* 

0-1 

0-1 

3-8 

Tracking efficiency 

0-1 and 11* 

0-1 and 9-11* 

9-11 and 1* 

9-11 and 0-1* 

8-14 and 1* 

Trigger efficiency 

0-1 and 2* 

0-1 and 2* 

2^ and 1* 

2^ and 0-1* 

2 and 1* 

Matching efficiency 

1* 

1* 

1 

1 

1 

pp stat. 

4+ 

5-12* 

6-21 

6-21 

7-11 

^J/w 

' ^ syst. 

8* 

7* 

5-6 and 6* 

5-6 and 6* 

5-6 and 6* 

^norm 

4+ 

4* 

4* 

4+ 

4+ 

(Taa) 

3-8 

3-6 

3* 

3-5* 

3* 

Centrality limits 

0-5 

0-3 

0 

0-2* 

0 

B.R. 

n/a 

n/a 

n/a 

1* 

n/a 


Table 2: Summary of the systematic uncertainties (in %) entering the J/y/ yield and/or /?aa calculation as a function of 
centrality, pj and y. Numbers with an asterisk correspond to the systematic uncertainties fully correlated as a function 
of the given variable. 



Fig. 5 : Differential yields of inclusive }/y in Pb-Pb collisions at = 2.76TeV as a function of pj for three 
centrality classes. Solid lines correspond to the results from the fit described in the text. 
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inclusive J/y yields with the function 


fipj)=Cx 


Pt 


i^ + iPj/PoV) 


(V) 


where C, po and n are free parameters. This function is commonly used to reproduce the J/t/r p^ distribution 
in hadronic collisions, see for instance |54-5^. Fit results for the three centrality classes are displayed as 
full lines in the figure. An excess over this function is revealed in the lowest pj interval (corresponding 
to 0 < Pt < 500MeV/c) for peripheral Pb-Pb collisions. It could be caused by a residual contribution 
from J/t/r coherent photo-production, which was measured in ulha-peripheral collisions 14211 . A quantitative 
measurement of this contribution in hadronic collisions is reported in 0]. Thus, in the most peripheral 
centrality class (40-90%) the fit was performed for pj > 500MeV/c and extrapolated down to zero (dotted 
line). In the 0-20% and 20-40% cenhality classes, no J/tp excess was observed and fits were performed 
down to zero pj. Asa cross-check, the same procedure as for the peripheral cenhality class was tested and 
the obtained results are fully compatible within uncertainties. 

Values of the mean hansverse momentum ((px)) and mean squared hansverse momentum ((px)) obtained 
from the fits are given in Tab.|3]as a function of centrality. The statistical (systematic) uncertainty is exhacted 
by fitting the px dishibution considering only the statistical (px-uncorrelated systematic) uncertainty of the 
measurement. For comparisons, the (px) and (p^) results from PHENIX were recomputed with the function 
defined by Eq.|7J adjusted in the measured px range and exhapolated to px = 8GeV/c to match our px range. 
These results are also given in Tab. [3] along with the measurement in pp collisions at ^/s = 2.76 TeV with 
updated uncertainties 071] . 


The (px) of inclusive J/t/r measured in pp and Pb-Pb collisions at = 2.76 TeV is shown in Fig. |6] (left 
side) as a function of (Apart). The error bars (open boxes) represent the statistical (systematic) uncertainties. 
A clear downward trend in (px) is observed when going from pp to the most central Pb-Pb collisions. The 
(px) decrease from peripheral (40-90%) to central (0-20%) collisions is significant, the two values being 
separated by more than 5a. These results are compared to the ones obtained by PHENIX in pp, Cu-Cu 
and Au-Au collisions at = 0.2 TeV. There is no evidence for a decreasing hend, contrary to what is 
observed in the AFICE measurement. 

In order to compare the evolution of (pp)a-a at different energies, one can form the variable taa defined as 


rAA = 


{P 


T/A—A 


{P' 


( 8 ) 


T/PP 


This variable was measured over the wide range of energies and colliding systems covered by NASO and 
PHENIX experiments. The comparison with the AFICE results is done in Fig. (right side). A very 
different (Apart) dependence is seen, especially when comparing Pb-Pb collisions at the SPS and the EHC. 
At the SPS energy of = 0.017 TeV the increase of the J/t/r (p^) with the cenhality of the collision 
was attributed to the Cronin effect 16IH . interpreted as an extra px kick due to multiple scatterings of the 
initial partons producing the J/t/r. At the EHC, a clear decrease of taa is observed as a function of (Apart). 
This behavior could be related to the onset of recombination phenomena and to the thermalization of charm 
quarks. Theoretical calculations 113, based on hansport models (described in the next section) are able 
to reproduce the taa at SPS, RHIC and EHC energies. They correlate the specific dependence of taa on 
collision cenhality with the increased importance of recombination effects in the J/t/r production mechanism 
at the EHC. 
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Px range 

y range 

Cenfralify 

(px) ± sfaf. ± syst. 

(Px) ± sfaf. ± sysf. 

(GeV/c) 



(GeV/c) 

(GeVVc^) 



Pb-Pb. 

y5;^ = 2.76TeV 


0-8 

2.5-4 

0-20% 

1.92 ± 0.02 ± 0.03 

5.17 ±0.12 ±0.16 

0-8 

2.5^ 

20-40% 

2.04 ± 0.02 ± 0.04 

5.83 ±0.11 ±0.17 

0.5-8 

2.5^ 

40-90% 

2.22 ± 0.03 ± 0.04 

6.72 ±0.14 ±0.20 



pp 

= 2.76TeV 


0-8 

2.5^ 

n/a 

2.28 ± 0.04 ± 0.03 

7.06 ±0.26 ±0.13 



PP 

= 0.2TeV [55] 


0-7 

1.2-2.2 

n/a 

1.61 ±0.01 ±0.012 

3.60 ± 0.06 ± 0.07 



Au-Au ^ 

t;];; = 0.2TeV [21J 


0-5 

1.2-2.2 

0-20% 

1.94 ±0.18 

5.79 ± 1.33 

0-6 

1.2-2.2 

20-40% 

1.87 ±0.07 

4.78 ± 0.34 

0-6 

1.2-2.2 

40-60% 

1.74 ±0.04 

4.19 ±0.27 

0-6 

1.2-2.2 

60-92% 

1.61 ±0.05 

3.87 ± 0.27 



Cu-Cu ^ 

^ = 0.2TeV|^] 


0-5 

1.2-2.2 

0-20% 

1.68 ±0.04 ±0.02 

3.79 ±0.25 ±0.11 

0-5 

1.2-2.2 

20-40% 

1.69 ±0.04 ±0.02 

3.71 ±0.18 ±0.08 

0-5 

1.2-2.2 

40-60% 

1.68 ±0.05 ±0.02 

3.91 ±0.30 ±0.11 

0-5 

1.2-2.2 

60-94% 

1.66 ±0.10 ±0.04 

4.13 ±0.64 ±0.24 


Table 3: Values of (px) and (p^) at various energies and colliding systems. The statistical and systematic uncertainties 
are quoted separately, except for PHENIX measurements in Au-Au collisions where the quadratic sum is given. If the 
measurement is not available or not used in the range 0 < px < 8 GeV/ c, the fit function is extrapolated down to 0 and 
up to 8GeV/c to compute (px) and (pj). 


9 Nuclear Modification Factor 


Some of the Raa results presented here were already published in ll27ll and are shown again in this section, 
where they are compared with model calculations and with results from previous experiments. They include 
the centrality dependence of Raa (Fig. ID, the px dependence of Raa for the full centrality range 0-90% 
and for the centrality class 0-20% (Fig. |9]top row) and the rapidity dependence of Raa (Fig.[l0|. The new 
results shown in this section include the centrality dependence of Raa for three pj intervals (Fig. [Hi and 
the px dependence of Raa for the centrality classes 20-40% and 40-90% (Fig. |9] bottom row). These new 
results were obtained using a slightly different trigger selection, as explained in section [3] The consistency 
of the results obtained with the two selections was verified. 


9.1 Centrality dependence of Raa 

Our measurement of the inclusive J/t/r Raa at = 2.76TeV in the range 2.5 < y < 4 and px < 8GeV/c 
is shown in Fig. [7] as a function of (Apart)- Statistical (uncorrelated systematic) uncertainties are represented 
by vertical error bars (open boxes). A global correlated systematic uncertainty affecting all the values by 
the same amount is quoted in the legend. The same convention is applied in the following figures, unless 
ofherwise specified. The J/y ^aa in the cenfralify class 0-90% (corresponding fo (Apart) ~ 124, see 
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Fig. 6: Mean transverse momentum (px) measured by ALICE 0371 ISTfl and PHENIX 02 iL 1551 l58ll as a function of 


the number of participant nucleons (left), taa measured by NA50 PHENIX and ALICE and compared to model 
calculations ES, as a function of the number of participant nucleons (right). 


Tab.lB is ~ 0.58 ± 0.01 (stat.) zb0.09(syst.), indicating a clear J/t/r suppression. This suppression is 


significantly less pronounced than that observed at lower energy in PHENIX in a similar kinematic range, 
as previously discussed in lE El- Eor (Apart) larger than 70, corresponding to the 50% most central Pb-Pb 
collisions, the J/ty /?aa is consistent with being constant, within uncertainties. Such behavior was not 
observed in heavy ion collisions at lower energies (SPS, RHIC), where R^a is continuously decreasing as a 
function of centrality. 


The impact of non-prompt J/t// on the inclusive Raa analysis was studied. The Raa of prompt J/ty is 
estimated (see Eq. O to be about 7% larger than the inclusive J/t// Raa if the beauty component is fully 
suppressed. In the other extreme case, where the B-meson production is not affected by the medium and 
scales with the number of binary collisions, i.e. 7?^™ prompt _ prompt J/ty would be about 

6% smaller in central collisions and about 1 % smaller in peripheral collisions. The excess in the inclusive 
J/Y yield observed at very low lliS also influences the Raa in the most peripheral collisions. A large 
fraction of this contribution (about 75% as explained in section lUl can be removed by selecting J/y with 
a Pt higher than 0.3GeV/c. Assuming that the hadronic 3/y Raa in the ranges 0 < px < 0.3GeV/c and 
0.3 < Pt < 8GeV/c are the same, it becomes possible to estimate the impact of the J/y photo-production 
on the inclusive Raa- In the centrality classes 60-70%, 70-80% and 80-90%, the hadronic J/y Raa would 
be about 5%, 11% and 25% lower, respectively. Extreme hypotheses were made to define upper and lower 
limits, represented with brackets on the Eigs. |7j [8]and|9] The upper limit calculation assumes no J/ty 
from photo-production thus the inclusive measurement only contains hadronic production. The lower limit 
assumes that i) all J/i/r produced with a px smaller than 0.3GeV/c originate from photo-production and ii) 
the efficiency of the 0.3GeV/c px selection is reduced from 75% to 60% (corresponding to an increase by a 
factor two of the J/y photo-production above 0.3GeV/c). 


The comparison with theoretical models, shown on the right-hand side of Eig. |7l helps in the interpretation 
of the large difference observed between the PHENIX and the AEICE results. 


The Statistical Hadronization Model (SHM) i62ll assumes deconfinement and thermal equilibration of the 
bulk of the cc pairs. Charmonium production occurs at the phase boundary via the statistical hadronization 
of charm quarks. The prediction is given for two values of the charm cross section dacc/dy = 0.15 and 0.25 
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Fig. 7: Inclusive J/y/ /?aa as a function of the number of participant nucleons measured in Pb-Pb collisions at 
= 2.76TeV compared to the PHENIX measurement in Au-Au collisions at = 0.2TeV ll^ (left) 
and to theoretical models |T3, S W], which all include a J/i// regeneration component (right). The brackets 
shown in the three most peripheral centrality classes on the right figure quantify the possible range of variation of the 
hadronic J /y Raa for two extreme hypotheses on the photo-production contamination in the inclusive measurement, 
see text for details. 

mb at forward rapidity. These values are derived from the measured charm cross section in pp collisions at 
y/s = 2.76 and 7TeV ITsIl bracketing the expectation for gluon shadowing in the Pb-nucleus between 0.6 
and 1.0. Production of non-prompt J/y from decays of B-mesons is not considered. 


The two transport models from Zhao (TMl) 111311 and Zhou (TM2) ll60ll mainly differ in the rate equa¬ 
tion controlling the J/y dissociation and regeneration. In TMl, shadowing is implemented via a simple 
parametrization, leading to a 30% suppression in the most central Pb-Pb collisions. The charm cross sec¬ 
tion is assumed to be dacc/dy ss 0.5 mb at forward rapidity, the fraction of J/y from beauty hadrons to be 
10% and no b-quenching is introduced in the calculation. This model is presented as a band connecting the 
results obtained with (lower limit) and without (upper limit) shadowing and is interpreted by the authors 
as the uncertainty of the prediction. In TM2, the shadowing is given by the EKS98 parametrization 16411 . 
The charm cross section is taken in the range dacc/dy 0.4 — 0.5 mb at forward rapidity; the calculations 
for these two values provide the lower and upper limits of the band displayed in the figure. The fraction of 
J/y from beauty hadrons is assumed to be 10% with a b-quenching of 0.8, increased to 0.4 for pj above 
5GeV/c. 


The Comover Interaction Model (CIM) ll63ll implements shadowing, interaction with a co-moving dense 
partonic medium and recombination effects. The shadowing is calculated within the Glauber-Gribov theory 
making use of the generalized Schwimmer model of multiple scattering. The J/y dissociation cross section 
due to comover interaction is taken as aco = 0.65 mb from low-energy data. Recombination effects are 
included by adding a gain term proportional to aco and to the number of c and c quarks, thus no additional 
parameter is added to the model. The charm cross section dacc/dy at forward rapidity is taken in the range 
0.4 to 0.6 mb, which gives respectively the lower and upper limits of the calculation. Production of non¬ 
prompt J/t// is not considered. 

To match our J/y Raa results, all models above need to include in their calculation a sizeable J/y production 
from deconfined c and c quarks. 
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Fig. 8: Inclusive J/i/ Raa as a function of the number of participant nucleons measured in Pb-Pb collisions at 
= 2.76TeV for three pj ranges (0-2, 2-5 and 5-8GeV/c) and comparisons of the lowest and highest px range 
to the transport and to the comover interaction models ifoi 1^ 1^ . The brackets quantify the possible range of 
variation of the hadronic J/ y/ Raa for two extreme hypotheses on the photo-production contamination in the inclusive 
measurement. 

A different test of these models was carried out by studying the J/yr Raa centrality dependence in px 
intervals. Figure |8] displays the measurement of the inclusive J/yr Raa as a function of the number of 
participant nucleons measured in Pb-Pb collisions at = 2.76TeV for the three pj ranges 0-2, 2-5 and 
5-8GeV/c. The uncorrelated systematic uncertainties shown at each point were separated into uncorrelated 
as a function of centrality (open boxes) and fully correlated as a function of centrality but uncorrelated as 
a function of pj (shaded areas). For (Apart) ^ 150, the low Pt^Iy ^aa is significantly larger than the mid 
and high px ones. In the most central bin, the Raa values corresponding to the lowest and the highest px 
are separated by 3.9a. For (Apart) ^ 150, the centrality dependence exhibits similar trends for the 2-5 and 
5-8GeV/c ranges, while the most peripheral ((Apart) ~ 20) Raa measurement in the low px (0-2 GeV/c) 
range appears to deviate from the others. However, the J/t/r yield excess observed at very low px may have 
a sizable effect in the 0-2 GeV/c interval. In the centrality classes 40-50%, 50-60% and 60-90%, based 
on the same assumptions made for the 0 < px < 8GeV/c case, the hadronic J/i/r Raa would be about 5%, 
6% and 18% lower, respectively. Due to the increase of the non-prompt J/t/r component at large px, the 
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difference between the measured inclusive J/t/r and the prompt J/t/r increases with pj. If the 
beauty contribution is fully (not) suppressed, /?aa of prompt J/t/r is estimated to be 6%, 8% and 11% larger 
(0-3%, 3-10% and 7-30% smaller, depending on centrality) for the pj ranges 0-2, 2-5 and 5-8GeV/c, 
respectively. 

Calculations from the transport models and the comover interaction model are plotted on top of the results 
shown in Fig. [S] For the most peripheral collisions ((Apart) ^ 100), the models cannot correctly reproduce 
the Raa centrality dependence for both the low and high pf ranges. For the most central collisions ((Apart) ^ 
100), the Raa centrality dependence for high PjJ/y is reasonably reproduced by all models. Concerning the 
low Pj range in the most central events, the measurement is compatible with the upper side of the theoretical 
uncertainty band from the CIM and TM2 models. For these models, it corresponds to the highest value for 
dacc/dy, 0.6 and 0.5 mb respectively. 


9.2 Transverse momentum dependence of Raa 

The Pj dependence of the inclusive J/t/r Raa in the rapidity range 2.5 < y < 4 is shown in Fig.|9]for the full 
centrality range 0-90% and for three centrality classes 0-20% |27], 20-40% and 40-90%. In Fig. |9] top left 
comer, the inclusive J/y/ Raa in the centrality class 0-90% shows a decrease of about 50% from low to high 
Pj. At low Pj, the measurement is close to 0.8 showing very little suppression. At high pj, our Raa value 
is similar to that of CMS Q. They measured, in the different rapidity range 1.6 < |y| < 2.4, an inclusive 
J/Y Raa = 0.41 ±0.05±0.04 for 3 < px < 30GeV/c. The corresponding mean pj is 6.27 GeV/c. When 
beauty contribution is fully (not) suppressed the prompt J/tp Raa is expected to be 5% larger (2% smaller) 
for Pj < IGeV/c and 17% larger (30% smaller) for 6 < px < 8GeV/c. 


The transport model calculations TMl 11311 and TM2 ihOT are also shown in Fig.|9] Both models reproduce 
reasonably well the 0-90% centrality measurement at high pj. At low px, TMl reproduces rather well our 
measurement, while the data points sit on the upper limit of the TM2 calculation. One can also appreciate the 
relative contributions of the primordial (from the initial hard parton scattering) and regenerated (from coales¬ 
cence of c and c quarks in the deconfined medium) components in these two calculations. The contribution 
of regenerated J/t/r is concentrated at low px and its relative fraction with respect to the initial production 
differs between the models. In TMl, it is of the same order of the primordial J/tp production, which is about 
constant over the full px range. In TM2, the regenerated J/tp contribution is almost three times larger than 
the primordial one in the lowest px interval. For px > 5GeV/c, only the primordial production remains. 

In the other panels of Fig. |9j the AFICE measurements are compared to those from PHENIX in Au-Au col¬ 


lisions at = 0.2TeV for the 0-20%, 20-40% and 40-60% centrality classes 12111 . For px < IGeV/c, 


for all centrality ranges, the prompt J/tp Raa is expected to be 5% larger (2% smaller) when the beauty 
contribution is fully (not) suppressed. For 6 < px < 8GeV/c the effect is much larger: if the beauty contri¬ 
bution is fully suppressed, the prompt J/t/r Raa would be 17% larger in all centrality ranges. If the beauty 
contribution is not suppressed, the prompt J/t/r Raa would be 44%, 15% and 8% smaller in the central¬ 
ity ranges 0-20%, 20-40% and 40-90%, respectively. The very low px excess in the inclusive J/tp yield 
mentioned before has a non-negligible impact in the 0 < px < IGeV/c range in the most peripheral cen¬ 
trality class 40-90%: following the same assumptions made for the 0 < px < 8GeV/c case, the estimated 
hadronic i/Y Raa would be about 22% lower. In the most central collisions (0-20%), the inclusive J/t/r 
Raa at low px is almost four times larger in Pb-Pb collisions at = 2.76TeV than in Au-Au collisions 
“ = 0.2 TeV. This difference cannot be explained only by the possible change in the size of the CNM 


at 


effects that can be expected due to the different rapidity coverage and collision energy between the two 
measurements. Such a behavior, on the other hand, is expected by all the recombination models described 
in the previous section. The same trend is observed in the centrality class 20-40%, where the large dif¬ 
ference between the PHENIX and AFICE results observed at low pj vanishes at high pj. Concerning the 
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Fig. 9: Inclusive J/i/ Raa as a function of the ]/y Pi for 2.5 < y < 4 in the centrality class 0-90% El compared to 
transport models |fH Ibol] (top left). The comparison is done with PHENIX results |21] and ttansport models in the 
0-20% 127] (top right), 20-40% (bottom left) and 40-90% (bottom right) centrality classes. The brackets shown in 
the lowest interval for the centrality class 40-90% quantify the possible range of variation of the hadronic }/y Raa 
for two extreme hypotheses on the photo-production contamination in the inclusive measurement. Upper limits from 
PHENIX at high p^ are not represented. 


most peripheral collisions, the inclusive J/y Raa is still slightly larger for ALICE results at low pj. How¬ 
ever, here the comparison between the two experiments is done with different centrality classes, 40-90% 
(ALICE) and 40-60% (PHENIX), so that a hrm conclusion, also because of the uncertainty size, cannot 
be drawn. Transport model calculations for Pb-Pb collisions at = 2.76TeV are also presented for the 
0-20%, 20-40% and 40-90% centrality classes. TMl shows a good agreement with the measurements in 
the 0-20% and 20-40% centrality classes, while TM2 tends to underestimate the data for pj < 5GeV/c. 
In the most peripheral centrality class (40-90%), the two models follow signihcantly different trends, but 
the uncertainties from the measurement are too large to discriminate them. However, if the very low pj 
excess is taken into account, a rather flat px dependence of the Raa is expected, pushing our measurement 
aside from TMl calculations in this specific range. Eor the high p^ region, both models reproduce well the 
experimental results in all the centrality classes. 
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9.3 Rapidity dependence of /?aa 


The rapidity dependence of the inclusive J/t/r /?aa in Pb-Pb collisions 11271] is shown in Fig.[T0l The inclusive 
J/Y Raa measured in the rapidity range |y| < 0.8 is about 0.7, consistent with the value measured at y ~ 3. 
From y ~ 3 to y ~ 4, the J/y Raa shows a decreasing trend leading to a drop of about 40%. The influence 
of non-prompt J/y/ on this result is small, as the prompt J/y Raa is expected to be only 8% larger (5% 
smaller) for 2.5 < y < 2.75 and 6% larger (9% smaller) for 3.75 < y < 4 if the beauty contribution is fully 
(not) suppressed. 



Fig. 10: Inclusive J/y/ Raa as a function of the J/y rapidity measured in Pb-Pb collisions at = 2.76 TeV ll27ll . 
compared to theoretical calculations of CNM effects due to shadowing and/or coherent energy loss 11651467 1. 


The Pb-Pb measurements are compared to theoretical calculations, which only consider shadowing and 
coherent energy loss. The break-up of the cc pair and nuclear absorption are not taken into account in any 
of the models. Shadowing only predictions are made within the Color Singlet Model at Feading Order |^] 
and the Color Evaporation Model at Next to Feading Order 1^ . with the EKS98 \M] and the EPS09 |^] 
parametrizations of the nPDF, respectively. For EKS98 (EPS09) the upper and lower limits correspond to 
the uncertainty in the factorization scale (uncertainty of the nPDF). Finally, a theoretical prediction, which 
includes a contribution from coherent parton energy loss processes in addition to EPS09 shadowing is 
also shown. All models show a fair agreement with our measurements over a wide rapidity range, |y| <3. 
If the amplitude of CNM effects is correctly given by the calculations shown in Fig.[T0l the observed J/y 
suppression due to CNM effects could be as large as 40%. Moreover, if an additional J/y suppression occurs 
in the hot nuclear matter (as expected from lower energy measurement and observed at high pj by both CMS 
and AFICE), other mechanisms compensating this suppression are needed to explain the Raa measurements. 
Figure [TO] supports this scenario, where suppression effects in hot matter are qualitatively counterbalanced 
by recombination. This is indeed what is expected from all models featuring recombination discussed in 
this paper. At higher rapidity, for |y| > 3, the models implementing only CNM effects tend to deviate from 
the data, although the one combining shadowing with coherent energy loss seems to match the decreasing 
trend of the Raa better. Such a decrease of the Raa values can also be explained by recombination models, 
where a reduction of the recombination effects is expected with increasing rapidity, due to the decrease of 
dacc/dy. 
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10 [y{2S)/J/y] ratio 


The ratio between inclusive Yi^S) and J/y yields measured in Pb-Pb collisions at = 2.76TeV is 
shown in the left side of Fig.[TT]as a function of (Apart)- In the interval pj < 3GeV/c, the ¥(2S) signal was 
extracted in three centrality classes (20-40%, 40-60% and 60-90%) while only the 95% confidence level 
upper limit was established for the centrality class 0-20%. At higher in the interval 3 < pj < 8GeV/c, 
the yield of ¥(2S) could not be extracted and the 95% confidence level upper limit is shown for the 0-20% 
and 20-60% most central collisions. 
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Fig. 11: Inclusive [v/(2S)/J/v/] ratio measured as a function of (Apart) in Pb-Pb collisions at = 2.76TeV for 
two pr intervals, compared to NA50 results lE^ and to a theoretical calculation iQ (left). Double ratio, as a function 
of (Apart), between the V^(2S) and J/y measured in Pb-Pb at ^^777 = 2.76TeV and pp collisions at = 7TeV, 
compared to theoretical calculations o (right). 


Our results are compared to the corresponding rneasurement at SPS energy (^/S(77 = 0.017 TeV), performed 
in a region close to mid-rapidity (0 <y < 1) 12911 . Within the rather large uncertainties of our measurement, 
no clear ^/S777 or y-dependence can be seen, in agreement with expectations from the SHM ifl^ . Prediction 
from the SHM for the prompt [¥(2S)/J/¥]pb-Pb '"^ho at yA777 = 2.76TeV in our rapidity domain is also 
reported in Fig. [TTI 


The double ratio [t/r(2S)/J/i/r]p|^_p|^/ [y{2S)/J/y]^^ is shown as a function of (Apart) in the right-hand side 
of Fig. [TT] Statistical uncertainties (including those coming from Pb-Pb and from the normalization to pp) 
are shown as vertical bars, while systematic uncertainties are shown as open boxes. The results do not 
allow a firm conclusion since statisfical fluctuations inside one standard deviation allow our data points to 
range between very low double ratios (strong Yi^-S) suppression with respect to J/y) to values higher than 
unity (less Yi^-S) suppression with respect to J/y)- Nevertheless, the limit set for the lowest p^ bin for the 
0-20% most central collisions points to a larger suppression of the ¥(2S) in that region. A transport model 
calculation tflTII for inclusive y{2S) and J/y production is shown for the two Pt intervals considered. The 
theoretical uncertainty band is due to different choices of the quenching factor for the b-quark. A qualitative 
agreement can be appreciated for both pj intervals. 


CMS has measured the double ratio [i/r(2S)/J/i/r]p(j_p^ / [i/r(2S)/J/i/r]pp dependence on centrality OOll for 
prompt ¥(2S) and J/y- In the rapidity and transverse momentum intervals 1.6 < |y| < 2.4 and 3 < px < 
30GeV/c and for the 0-20% most central collisions, a double ratio of 2.31 ±0.53(stat.) ±0.37(stat.) ± 
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0.15(pp) is obtained. This result sits at the upper edge of our confidence limit estimated in the same cen¬ 
trality range for 2.5 <y <4 and 3 < Pt < 8GeV/c. In more peripheral collisions, CMS results fall inside 
the limits given by this analysis. 


It is worth underlying that our result is for inclusive J/t/r production, while SHM predictions and 

CMS results are for prompt charmonia production. The impact of the B-mesons feed-down on the ratio was 
extensively studied in ilTII . showing a very strong influence of the non-prompt V^(2S) component on the 
final result. According to this study, removing this non-prompt contribution would lead to a significantly 
lower double ratio at high (A^pait): in the 0 < px < 3GeV/c bin a 60% decrease is expected, while in the 
3 < Pt < 8GeV/c bin the effect could be even stronger, leading to a 80% decrease. 


11 Conclusions 


We have presented a study of J/tp and t//^(2S) production in Pb-Pb collisions at = 2.76TeV in the 
transverse momentum and rapidity ranges px < 8GeV/c and 2.5 < y < 4. This analysis was carried out in 
the muon spectrometer system, whose tracking and triggering capabilities were described in detail. 

The [t/r(2S)/J/t/r]p[,_p|, ratio was measured in two px ranges as a function of centrality. In some intervals, 
only the 95% confidence level upper limits could be obtained. The suppression pattern of the V^(2S) is 
compatible with that of the J/y in most of the centrality and px intervals studied. The large uncertainties 
leave open the possibility of strong enhancement or suppression factors. An accurate t/r(2S) measurement 
in Pb-Pb would require significantly more statistics than the one presented in this analysis. 

The J/y signal was extracted as a function of px, y and the collision centrality. We have computed the 
J/Y (Pt) and (px)- The J/y (Pt) in Pb-Pb collisions decreases significantly (5a effect) from peripheral to 
central collisions. In addition we have studied taa defined as the ratio of the J/y (Px) nieasured in Pb-Pb 
and pp collisions at the same energy. The taa exhibits a clear decrease as a function of centrality for Pb-Pb 
collisions. 


The nuclear modification factor, Raa, of inclusive J/y was measured as a function of centrality. A constant 
suppression of about 40% was observed for (Apart) larger than 70 12711 . New studies of the J/y suppression 
pattern as a function of centrality for three px ranges were presented. Above (Apart) ~ 150, the low px J/y 
Raa clearly differs from the high PjJ/y ^aa and is about three times larger for (Apart) > 250, corresponding 
to a 3.9a separation. Complementary to this, the px dependence of the suppression pattern was analysed 
for the different centrality classes. An increase of the inclusive J/y Raa with decreasing px is observed 
below 5 GeV/c in the most central Pb-Pb collisions (0-20%), while no significant px dependence is seen 
in the most peripheral collisions (40-90%). As a function of rapidity, the results published in 12711 show 
compatible Raa values for |y| ^ 0.9 and 2.5 <y< 3. For larger rapidity, a decreasing trend is visible. 


Comparisons of the vaa and Raa measured in ALICE with lower energy experiments show significant 
differences. The decreasing trend of taa observed as a function of centrality is opposite to NA50 and 
PHENIX measurements. The Raa in the most central collisions is three times larger than the one measured 
by PHENIX, and the difference reaches a factor four in the px region below IGeV/c. If the suppression 
sources observed at lower energies, which were related to color screening in hot nuclear matter on top of 
CNM effects, are still present at the LHC, then other mechanisms compensating the J/y suppression are 
needed to explain the ALICE measurements. This conclusion is further substantiated, in the region |y| < 3, 
by the comparison of the inclusive J/y Raa measurements as a function of y to models implementing only 
CNM effects, which shows a qualitative agreement. 

The inclusive J/y ^'aa and Raa measurements were also compared to various theoretical calculations in¬ 
cluding hot and cold nuclear matter effects. The hadronic part of the J/y Raa was estimated when needed 
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to allow for a direct comparison to models, which do not implement the J/t/r production mechanism at the 
origin of the observed very low excess 11401] . All these models feature a full or partial J/t/r production 
from charm quarks recombination and are in fair agreement with the experimental results. The transport 
models considered in this paper are also able to generate an amount of J/t/r elliptic flow comparable to the 
one measured in ALICE 16911 . The double differential studies of the inclusive J/yr /?aa as a function of 
centrality and pj brings new constraints to the models. Reproducing the suppression pattern in peripheral 
collisions for both low and high px is challenging for all models. Some tensions also appear in describ¬ 
ing the /?aa evolution at low pj for all centrality classes. However, the uncertainties on the measurements 
on one side, and on the CNM and dace/dj in the theoretical calculations on the other side, do not allow for 
drawing a firm conclusion. The large uncertainties on the model predictions also show the limit of using 
the Raa ns an observable to measure the J/t/r suppression due to hot medium effects. Ideally one should, in 
Pb-Pb collisions, compare the J/y/ production to the charm production to cancel out the cold nuclear matter 
effects affecting the initial cc dynamics. However, the measurement of the charm cross section in Pb-Pb 
collisions is very ambitious and still remains to be done at the LHC. 


To summarize, the J/y/ taa and Raa measured in Pb-Pb collisions at = 2.76TeV show a new behavior 

with respect to measurements made at lower energies. In addition to the strong J/y/ suppression observed 
at high Pj, ALICE results show that at low px a new contribution is necessary to explain the data. In all 
available model calculations, this contribution is related to a recombination mechanism of charm quarks. 
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A Data tables 


This appendix provides all the numerical values obtained in this analysis. 


The inclusive J/t/r differential pj yields in Pb-Pb in centrality classes are given in Tab. lA.ll Tables [A2l to 
IA.5[ present the inclusive J/yr /?aa and associated Pb-Pb yields as a function of centrality for 2.5 <y <4.0 
and four pj ranges, p^ < 8GeV/c, pj < 2GeV/c, 2 < pj < 5GeV/c and 5 < px < 8GeV/c. Tables IAT^ 
to lA. 101 show the px dependence of the inclusive J/y Raa and associated Pb-Pb yields for the centrality 
classes 0-20%, 20-40%, 0-40%, 40-90% and 0-90%. Table I A. lH shows the y dependence of the inclusive 
J/Y Raa and associated Pb-Pb yields for the centrality class 0-90% in the px range px < 8GeV/c. Then, 
the inclusive J/y Raa results with a low px cut at 0.3 GeV/c are presented. The reference pp cross section 
needed to build the Raa was extracted with the method described in i40n . The inclusive J/y Raa centrality 
dependence for 2.5 < y < 4 in the px ranges 0.3 < px < 8GeV/c and 0.3 < px < 2GeV/c is shown in 
Tab. lA. 121 The inclusive J/y Raa in the px range 0.3 < px < IGeV/c for 2.5 < y < 4 in four centrality 
classes 0-90%, 0-20%, 20-40% and 40-90% is given in Tab. IA.131 Einally, Tab. |AT4] presents the inclu¬ 
sive [tp(2S)/J/v/']pb_pb [V^(2S)/J/t/^]pb_pb / [t/r(2S)/J/tp]pp ratios as a function of centrality for the px 

intervals px < 3GeV/c and 3 < px < 8GeV/c. 
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d%^/dydpj{GtY/c)-^ X 103 


Pt (GeV/c) 

0-20% 

20-40% 

40-90% 

0.0-0.5 

3.253 ±0.386 ±0.446 

1.366±0.081 ±0.165 

0.257 ±0.017 ±0.031 

0.5-1.0 

8.012 ±0.487 ±1.087 

2.571 ±0.199±0.310 

0.346 ± 0.024 ±0.042 

1.0-1.5 

9.909 ±0.603 ±1.149 

3.494 ±0.255 ±0.388 

0.533 ±0.030 ±0.061 

1.5-2.0 

8.193 ±0.505 ±0.907 

2.907 ±0.194 ±0.320 

0.493 ±0.037 ±0.053 

2.0-2.5 

6.342 ±0.401 ±0.701 

2.371 ±0.164 ±0.260 

0.441 ±0.034 ±0.049 

2.5-3.0 

4.759 ±0.316 ±0.542 

1.997 ±0.134 ±0.227 

0.270 ±0.020 ±0.029 

3.0-3.5 

2.735 ±0.183 ±0.290 

1.313±0.087±0.151 

0.222 ±0.016 ±0.023 

3.5-4.0 

1.876±0.134±0.201 

0.874 ±0.068 ±0.092 

0.174±0.013±0.018 

4.0-4.5 

1.075 ±0.098 ±0.109 

0.483 ±0.037 ±0.048 

0.108 ± 0.009 ±0.011 

4.5-5.0 

0.731 ±0.069 ±0.073 

0.339 ±0.030 ±0.033 

0.076 ± 0.007 ±0.007 

5.0-5.5 

0.453 ±0.047 ±0.045 

0.263 ±0.023 ±0.026 

0.042 ± 0.005 ±0.004 

5.5-6.0 

0.345 ±0.039 ±0.046 

0.132±0.016±0.014 

0.028 ± 0.004 ±0.003 

6.0-8.0 

0.099 ±0.009 ±0.010 

0.068 ± 0.005 ±0.007 

0.012 ±0.001 ±0.001 


Table A.l: Inclusive J/y/ yields (as defined by Eq. [T]i in pj intervals for the 0-20%, 20^0% and 40-90% most 
central Pb-Pb collisions. The rapidity range is 2.5<y<4. Statistical and systematic uncertainties are also reported as 
d^Fj/^/dydpT ± statistical uncertainty ± systematic uncertainty. A global systematic uncertainty of 4% affects all the 
values. A 2%, 1% and 2% systematic uncertainty, independent of p^, affects the centrality classes 0-20%, 20^0% 
and 40-90%, respectively. 


Centrality 

7?aa± (stat.) ± (syst.) [27] 

Ej/^± (stat.) ± (syst.) x 10^ 

0-10% 

0.557 ±0.019 ±0.024 

43.095 ±1.454 ±1.049 

10-20% 

0.573 ±0.020 ±0.022 

27.212 ±0.979 ±0.501 

20-30% 

0.598 ±0.022 ±0.020 

17.409±0.638±0.188 

30-40% 

0.577 ± 0.024 ±0.025 

9.671 ±0.406 ±0.211 

40-50% 

0.609 ±0.028 ±0.030 

5.413 ± 0.247 ±0.041 

50-60% 

0.725 ±0.036 ±0.043 

3.246 ±0.160 ±0.050 

60-70% 

0.839 ±0.041 ±0.058 

1.677 ±0.083 ±0.024 

70-80% 

0.849 ±0.063 ±0.068 

0.701 ±0.051 ±0.014 

80-90% 

1.094 ±0.106 ±0.104 

0.362 ±0.033 ±0.008 


Table A.l: Inclusive J/y Raa and Pb-Pb yields as a function of centrality, for pi < 8GeV/c and 2.5 < y < 4.0. 
Statistical and systematic uncertainties are also reported. A global systematic uncertainty of 15% (12%) affects all the 
Raa (yields) values. 
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Centrality 

/?AA±(stat.)±(syst.) 

}j/^± (stat.) ± (syst.) X 10^ 

0-10% 

0.732 ±0.034 ±0.041 

27.932 ±1.302 ±1.282 

10-20% 

0.733 ±0.035 ±0.028 

17.159±0.824±0.383 

20-30% 

0.715 ±0.038 ±0.024 

10.113±0.541 ±0.115 

30-40% 

0.678 ±0.040 ±0.033 

5.516±0.322±0.182 

40-50% 

0.641 ±0.044 ±0.032 

2.789 ±0.190 ±0.064 

50-60% 

0.839 ±0.048 ±0.056 

1.799 ±0.103 ±0.070 

60-90% 

1.104 ± 0.064 ±0.078 

0.559 ±0.032 ±0.016 


Table A.3: Inclusive J/y/ Raa and Pb—Pb yields as a function of centrality, for pj < 2GeV/c and 2.5 < y < 4.0. 
Statistical and systematic uncertainties are also reported. A global systematic uncertainty of 15% (12%) affects all the 
Raa (yields) values. 


Centrality 

/?AA±(stat.)±(syst.) 

}j/^± (stat.) ± (syst.) X 10^ 

0-10% 

0.425 ±0.019 ±0.017 

15.540±0.681 ±0.379 

10-20% 

0.461 ±0.019 ±0.016 

10.336±0.431 ±0.168 

20-30% 

0.529 ±0.022 ±0.018 

7.164 ±0.293 ±0.106 

30^0% 

0.498 ±0.025 ±0.027 

3.879±0.194±0.153 

40-50% 

0.595 ±0.030 ±0.029 

2.481 ±0.126 ±0.049 

50-60% 

0.675 ±0.042 ±0.041 

1.386 ±0.085 ±0.037 

60-90% 

0.722 ±0.044 ±0.050 

0.350 ±0.021 ±0.009 


Table A.4: Inclusive J/y/ Raa and Pb-Pb yields as a function of centrality, for 2 < pj < 5GeV/c and 2.5 < y < 4.0. 
Statistical and systematic uncertainties are also reported. A global systematic uncertainty of 14% (11%) affects all the 
Raa (yields) values. 


Centrality 

/?AA±(stat.)±(syst.) 

}j/^± (stat.) ± (syst.) X 10^ 

0-10% 

0.280 ±0.021 ±0.011 

1.093 ±0.081 ±0.027 

10-20% 

0.282 ±0.027 ±0.011 

0.677 ±0.064 ±0.016 

20-30% 

0.410 ±0.029 ±0.013 

0.594 ±0.042 ±0.006 

30^0% 

0.540 ±0.039 ±0.024 

0.449 ±0.033 ±0.012 

40-50% 

0.529 ±0.053 ±0.031 

0.236 ±0.024 ±0.009 

50-60% 

0.587 ±0.073 ±0.036 

0.129 ±0.016 ±0.004 

60-90% 

0.644 ±0.083 ±0.047 

0.033 ± 0.004 ±0.001 


Table A.5: Inclusive J/y/ Raa and Pb-Pb yields as a function of centrality, for 5 < pj < 8GeV/c and 2.5 < y < 4.0. 
Statistical and systematic uncertainties are also reported. A global systematic uncertainty of 18% (10%) affects all the 
Raa (yields) values. 
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Pt (GeV/c) 

/?AA±(stat.)±(syst.) [27] 

d^Fj/^/dydpx ± (stat.) ± (syst.)(GeV/c) ' x 10^ 

0-1 

0.803 ±0.084 ±0.113 

5.771 ±0.345 ±0.748 

1-2 

0.690 ±0.052 ±0.084 

9.134±0.411±0.987 

2-3 

0.505 ± 0.042 ±0.062 

5.539 ±0.284 ±0.604 

3^ 

0.381 ±0.037 ±0.046 

2.305 ±0.116 ±0.247 

4-5 

0.355 ±0.052 ±0.041 

0.905 ±0.068 ±0.090 

5-6 

0.282 ±0.048 ±0.032 

0.388 ±0.030 ±0.038 

6-8 

0.279 ±0.064 ±0.032 

0.100 ± 0.009 ±0.010 


Table A.6: Inclusive J /y Raa and Pb-Pb yields as a function of pi for the 0-20% centrality class and 2.5 <y<4.0. 
Statistical and systematic uncertainties are also reported. A global systematic uncertainty of 8% (4%) affects all the 
Raa (yields) values. 


Pt (GeV/c) 

/?AA±(stat.)±(syst.) 

d^Fj/^/dydpx ± (stat.) ± (syst.)(GeV/c) ^ x 10^ 

0-1 

0.733 ±0.080 ±0.097 

1.909 ±0.128 ±0.229 

1-2 

0.660 ±0.051 ±0.080 

3.189±0.154±0.344 

2-3 

0.543 ± 0.044 ±0.067 

2.167±0.106±0.238 

3^ 

0.493 ±0.048 ±0.060 

1.084±0.055±0.117 

4-5 

0.444 ±0.063 ±0.051 

0.411 ±0.027 ±0.040 

5-6 

0.399 ±0.067 ±0.045 

0.200 ±0.014 ±0.020 

6-8 

0.523 ±0.116 ±0.059 

0.068 ±0.005 ±0.007 


Table A.7: Inclusive ]/y Raa and Pb-Pb yields as a function of pi for the 20^0% centrality class and 2.5<y<4.0. 
Statistical and systematic uncertainties are also reported. A global systematic uncertainty of 8% (4%) affects all the 
Raa (yields) values. 


Pt (GeV/c) 

/?AA±(stat.)±(syst.) 

d^Fj/^/dydpx ± (stat.) ± (syst.)(GeV/c) ^ x 10^ 

0-1 

0.767 ± 0.074 ±0.105 

3.754 ±0.163 ±0.472 

1-2 

0.672 ± 0.046 ±0.082 

6.103 ±0.212 ±0.662 

2-3 

0.515 ±0.038 ±0.064 

3.865 ±0.134 ±0.428 

3^ 

0.411 ±0.038 ±0.049 

1.698±0.063±0.178 

4-5 

0.376 ±0.051 ±0.043 

0.655 ±0.033 ±0.064 

5-6 

0.315 ± 0.050 ±0.036 

0.296 ±0.016 ±0.029 

6-8 

0.340 ±0.075 ±0.038 

0.083 ±0.005 ±0.008 


Table A.8: Inclusive J /y Raa and Pb-Pb yields as a function of p^ for the 0^0% centrality class and 2.5 <y<4.0. 
Statistical and systematic uncertainties are also reported. A global systematic uncertainty of 8% (4%) affects all the 
Raa (yields) values. 
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Pt (GeV/c) 

/?AA±(stat.)±(syst.) 

d^Fj/^/dydj)T ± (stat.) ± (syst.)(GeV/c) ^ x 10^ 

0-1 

0.815±0.081 ±0.107 

0.305 ±0.015 ±0.036 

1-2 

0.732 ±0.059 ±0.090 

0.508 ±0.028 ±0.055 

2-3 

0.617 ±0.053 ±0.076 

0.354 ±0.020 ±0.038 

3^ 

0.627 ±0.062 ±0.074 

0.198±0.010±0.020 

4-5 

0.693 ±0.097 ±0.079 

0.092 ± 0.006 ±0.009 

5-6 

0.489 ±0.087 ±0.055 

0.035 ±0.003 ±0.003 

6-8 

0.646 ±0.150 ±0.072 

0.012 ±0.001 ±0.001 


Table A.9: Inclusive ]/y Raa and Pb-Pb yields as a function of for the 40-90% centrality class and 2.5<y<4.0. 
Statistical and systematic uncertainties are also reported. A global systematic uncertainty of 9% (4%) affects all the 
Raa (yields) values. 


Pt (GeV/c) 

Raa ± (stat.) ± (syst.) [27] 

d^Fj/^/dydj)T ± (stat.) ± (syst.)(GeV/c) ^ x 10^ 

0-1 

0.779 ±0.076 ±0.106 

1.857 ±0.081 ±0.230 

1-2 

0.677 ±0.047 ±0.083 

2.993 ±0.104 ±0.323 

2-3 

0.519 ±0.038 ±0.064 

1.896 ±0.064 ±0.206 

3^ 

0.425 ±0.039 ±0.051 

0.855 ±0.029 ±0.089 

4-5 

0.405 ± 0.054 ±0.047 

0.343 ±0.015 ±0.033 

5-6 

0.322 ±0.052 ±0.036 

0.147 ± 0.007 ±0.015 

6-8 

0.364 ±0.079 ±0.041 

0.043 ±0.002 ±0.004 


Table A.IO: Inclusive J/i/ Raa and Pb-Pb yields as a function of pi for the 0-90% centrality class and 2.5 <y <4.0. 
Statistical and systematic uncertainties are also reported. A global systematic uncertainty of 8% (4%) affects all the 
Raa (yields) values. 


F 

/?aa± (stat.) ± (syst.) [27] 

d^Fj/^/dyd;?T ± (stat.) ± (syst.)(GeV/c) ' x 10^ 

2.50-2.75 

0.631 ±0.087 ±0.088 

1.509±0.114±0.191 

2.75-3.00 

0.747 ±0.068 ±0.097 

1.387±0.058±0.162 

3.00-3.25 

0.632 ±0.048 ±0.094 

1.120±0.039±0.154 

3.25-3.50 

0.566 ± 0.044 ±0.088 

0.891 ±0.032 ±0.130 

3.50-3.75 

0.467 ±0.041 ±0.070 

0.733 ±0.025 ±0.101 

3.75^.00 

0.395 ±0.050 ±0.050 

0.528 ±0.029 ±0.058 


Table A.ll: Inclusive J/y/ Raa and Pb-Pb yields as a function of y for the 0-90% centrality class and Pt < 8GeV/c. 
Statistical and systematic uncertainties are also reported. A global systematic uncertainty of 8% (4%) affects all the 
Raa (yields) values. 
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/?AA±(stat.)±(syst.) 

Centrality 

0.3 <pi< 8GeV/c 

0.3 < pi<2G&\/c 

0-10% 

0.545 ±0.017 ±0.026 

0.745 ±0.041 ±0.042 

10-20% 

0.560 ±0.018 ±0.021 

0.736 ±0.036 ±0.028 

20-30% 

0.594 ± 0.020 ±0.020 

0.716 ±0.038 ±0.025 

30-40% 

0.570 ±0.021 ±0.025 

0.671 ±0.040 ±0.032 

40-50% 

0.592 ±0.025 ±0.029 

0.619 ±0.045 ±0.032 

50-60% 

0.715 ±0.033 ±0.044 

0.801 ±0.049 ±0.054 

60-70% 

0.805 ±0.043 ±0.057 

'I 

70-80% 

0.778 ±0.062 ±0.064 

^0.959 ±0.057 ±0.067 

80-90% 

0.887 ±0.097 ±0.088 

J 


Table A.12: Inclusive J/y/ /?aa as a function of centrality, for 0.3 < pi < 8GeV/c and 0.3 < p-r < 2GeV/c in the 
rapidity range 2.5<y<4.0. Statistical and systematic uncertainties are also reported. A global systematic uncertainty 
of 15% affects all the /?aa values. 


Centrality 

/?aa± (stat.) ± (syst.) for 0.3 < pi < IGeV/c 

0-90% 

0.775 ±0.057 ±0.113 

0-20% 

0.803 ±0.066 ±0.123 

20-40% 

0.733 ±0.067 ±0.103 

40-90% 

0.688 ±0.057 ±0.098 


Table A.13: Inclusive J/y/ /?aa for 2.5 < y < 4.0 in the centrality classes 0-90%, 0-20%, 20-40% and 40-90% for 
the lowest pi range when the 0.3GeV/c pi cut is applied. Statistical and systematic uncertainties are also reported. A 
global systematic uncertainty of 8%, 8%, 8% and 9% affect the I?aa values, respectively. 


Pi (GeV/c) 

Centrality 

[’/^(2S)/J/t//']pb_pb 

[t/r(2S)/J/t/r]pj^p, / [Yi2S)/J/Y]„ 

0-3 

0-20% 

< 0.012 (95% CL) 

< 0.65 (95% CL) 

0-3 

20-40% 

0.017 ±0.010 ±0.004 

0.86 ±0.51 ±0.23 

0-3 

40-60% 

0.013 ±0.012 ±0.006 

0.65 ±0.65 ±0.30 

0-3 

60-90% 

0.029 ±0.012 ±0.004 

1.49 ±0.62 ±0.27 

3-8 

0-20% 

< 0.046 (95% CL) 

< 1.71 (95% CL) 

3-8 

20-60% 

< 0.033 (95% CL) 

< 1.24 (95% CL) 


Table A.14: Inclusive [v/(2S)/J/vr]pb_pi, and [v/(2S)/J/vr]pb_pb / [’/(2S)/J/v^]pp ratios as a function of centrality for 
two Pi intervals. Statistical and systematic uncertainties are reported when the value is not given as an upper limit. 
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